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The bias-dependent resistance R{V) of NS-junctions is cal- 
culated using the Keldysh formalism in all orders of the trans- 
fer matrix element. We present a compact and simple for- 
mula for the Andreev current, that results from the coupling 
of electrons and holes on the normal side via the anomalous 
Green's function on the superconducting side. Using sim- 
ple BCS Nambu-Green's functions the well known Blonder- 
Tinkam-Klapwijk theory can be recovered. Incorporating the 
energy-dependent quasi-particle lifetime of the heavy fermions 
strongly reduces the Andreev-reflection signal. 



I. INTRODUCTION 

Point-contact spectroscopy (PCS) has been used to 
study the superconducting (SC) properties of heavy- 
fermion (HF) compounds, see e.g. Refs. flj-^l- The 
symmetry of the SC order parameter in HF compounds 
is still not known, and it was hoped that PCS is a useful 
tool to help clarify this question. 

UsuaUy, the SC anomalies of point contacts between 
HFSC and normal metals have been interpreted in terms 
of Andreev reflection (AR). But the observed spectra 
deviate significantly from that predicted by the simple 
quasi-classical BTK model 0. For example, the dou- 
bling of the electrical current due to AR at a bias voltage 
\eV\ = A is missing. In recent experiments on UPta 0] 
and URu2Si2 Q (that are good candidates for being per- 
formed in the ballistic limit as claimed by the authors) 
the size of the SC anomaly amounts to a small fraction of 
the total signal. We present here a microscopic approach 
using the Keldysh non-equilibrium method. It incorpo- 
rates the strongly energy-dependent lifetime in HF com- 
pounds, that is also responsible for the large A coefficient 
of the electrical resistivity at temperatures T <^T*, the 
lattice Kondo temperature. For vanishing self-energy the 
BTK result is recovered. 



II. THEORY 

We use a simplified model to describe a NS-junction. 
Let us assume the left side of the junction is described by 
the Hamiltonian Hl, the right by Hn, and the solution 
of both is known in terms of the one-particle equilibrium 
Green's function. The two leads are coupled by a tunnel 
Hamiltonian |^,^ 
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<L,R> 



h.c. 
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The applied bias voltage translates into a step-like chem- 
ical potential at the junction with the potential difference 
eV . Hereby the assumption of a ballistic point contact 
enters the theory. For a diffusive or thermal junction, 
the spatial dependence of the chemical potential has to 
be determined self-consistently. The total current 



^M = |lETr[r3(i£+,(r,r)-|t^y,,. 



L.R 



(2) 



is in general time-dependent. We introduce the Keldysh 
Nambu 2x2 occupational Green's function 



£L^(r,r')=* 



with 



« c[|(T)cflT(T') » « c{^{t) 4^(r') » 
< clx(t) CflT(T') > < CLi(r) 4^(r') > 
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In this Nambu representation the hopping matrix ele- 
ment picks up a time-dependent phase difference (/'(r) = 
(j)Q + tujq/2 between the two leads, with luq = 2eV/h. 

The current can be decomposed in a Fourier series 
of multiples of the fundamental frequency ujq'- I{t) = 
/„ exp(iT7ia;o) necessary for the proper description of 
a AC Josephson current in SS-junction. In a NS-junction, 
however, all but the /q component vanish since only |Tp 
enters in this case. Applying the Keldysh formalism de- 
veloped for the tunneling theory j^] yields the equations 
of motion in frequency space 



G;, (a>) = / (c.)tt {lu) = GUu) tt gr (^) 



(5) 



which can be solved, omitting the frequency argument 
and using 
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where the two hopping matrix elements have been ab- 
sorbed into the renormalized normal-state propagator 
|T|2 



g^(ti>) on the right side (R). The left 



side (L) becomes SC. The denominator indicates the re- 
summation of tunneling processes in infinite order. 
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The occupational-components are more complicated: 
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Here g denotes the equilibrium advanced (a), retarded 
(r), and occupational (±) Green's functions in absence 
of Ht, G the fully renormalized non-equilibrium Green's 
functions. Using Eqs. (||)-(||), the notation eito) for 
an electron, h{Lo) for a hole, and /(w) for the anoma- 
lous Nambu Green's function component, the total DC 
current /(V,t ~ 0) of Eq.(|^) decomposes into a quasi- 
particle current 



Iqp = 
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+25fte[(l - T'h'khl)T'hyi]{e+f^ 
and the Andreev current 



^Rfi)] (10) 
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T'e^R^l) 



with \Det[l-g^^ = 1(1 - T^h-j,hl)[l 

T^h^ReWl??- 

Iqp describes three different processes: an electron on 

the normal side R couples to an electron, a hole, or to the 
anomalous GF. The Andreev current couples an electron 
and a hole on the normal side under the influence of a 
Cooper pair. Its leading order is since it involves two 
quasi-particles crossing the junction. Note that since we 
resummed all orders, the hopping T can be tuned contin- 
uously from being much smaller than the bandwidth W , 
the tunneling regime, to 7 = ttT Np w 1, the metallic 
PCS regime. 

^<L R> denotes two solvable quasi-ld cases: 
1. Quantum point contact. The orifice involves only a 
few lattice sites. Translational invariance is broken and 
no momentum component is conserved at the boundary. 
It can be verified that using BCS propagators and identi- 
fying = ((1 — 7'^27)^ the BTK theory is reproduced 
by Eqs. (|l^) and (|l|). where J2<l r> counts the num- 
ber of channels. For example, the Andreev current per 
channel reads 
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(1 + 7^)2(1 - (t^/A)2) + 474(lj/A)2 
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(l+7^)VW~l + VI^/A|y 
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(13) 
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Obviously the transmission coefficient T(aj) is continuous 
at uj"^ = and its value ^ is independent of the effec- 
tive coupling 7. 

2. Two connected half spaces. In this case the momentum 
parallel to the fc||-plane is conserved. e(tj, fc||), /(w, fcy) 
and /(w,^) become /c^-summed Green's functions enter- 
ing Eq. ( p^ and (^ij) . The total current is the sum of all 
independent contributions. For example, in Eq. |ll| an 
electron is coupled to a hole of the same fc|| which trans- 
lates into the standard one-particle picture of A: ^ — fc. 
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FIG. 1. Resistance per conductance channel dV/dl vs bias 
voltage V as function of the dimensionless hopping parameter 
7. Inset: blow-up of the 7 = 0.65 curve. 



III. APPLICATION 

A microscopic theory of HFSC including quasi-particle 
lifetime effects is missing up to now. We have used in- 
stead the normal state Green's function for the Anderson 
lattice 



5(^) = 



1 



-E,(z) , -$ 
-^-S*(z) 



(15) 



with an isotropic order parameter $ = O.l-T* as a simple, 
not self-consistent approximation for the SC HF phase 

where 9ml],(u; - i5) = ^ ((f?)' + (^)'). This 

formalism can be easily extended by an anisotropic or- 
der parameter <I>fe(iLj„) using case 2. Only the conduction 
electrons contribute to the charge transport since the hy- 
bridisation is assumed /c-independent g . 

The results of Eqs. ^ and (|ll]) are shown in Fig. |^. 
Ignoring the life-time effects and the hybridisation in 
( p^ ) indeed reproduces the BTK result for 7 = 1. But for 
the HE case a strongly reduced 'V'-shape SC anomaly is 
found for an ideal metallic point contact (7 = 1). The 
'worse' the contact (7 < 1) the more BTK-like is the 
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shape of an 'ideal' contact, but its size is strongly reduced 
at zero bias. Because we have chosen a fixed order param- 
eter $ the renormalization due to the strongly reduced 
quasi-particle spectral weight leads to the 'shrinking' of 
the gap. The inset shows a blow-up of the 7 = 0.65 curve. 
This would come close to some of the spectra observed for 
UPta. The asymmetry reflects the asymmetry between 
the normal electrode (treated by a constant density of 
states) and the HFSC. 

IV. DISCUSSION 

According to our above calculations the almost flat 
spectra with small double-minimum structure are acces- 
sible only within a very narrow parameter range when 
couphng is reduced to about 7 « 0.65. This could be one 
of the reasons why those spectra are very hard to detect 
Note that even with an isotropic order parameter 
spectra different from the BTK predictions can be ob- 
tained using a energy dependent quasi-particle life-time. 

We ought to mention some drawbacks to our ap- 
proach: The energy-dependence of the quasi-particle 
lifetime requires a self-consistent determination of the 
spatial-dependent potential to model contacts in the eas- 
ily accessible resistance range (O.Sili? < 3^1). In an 
Eliashberg-type of calculation of the SC state the onset 
of quasi-particle lifetime effects is expected to be shifted 
self-consistently towards the gap-edges. And we cannot 
explain the often found enhancement of the ballistic zero- 
bias resistance. 

A recent systematic study on the size of the SC anoma- 
lies clearly shows a scaling with the orifice radius rather 
than with the area This could indicate that a NNS- 
junction is seen in most of the experiments with a NS- 
boundary inside the HF material. Such a scenario would 
be also favoured in the case of an anisotropic order pa- 
rameter, since the surface then acts as a pair-braking 
boundary. A reminiscence of SC is found in the normal- 
state HF part of the junction, since it is known - or can 
be easily derived using Eq.(^ - that the non-equilibrium 
renormalized normal side has a gap at lu = in its spec- 
tral function induced by the superconductor it is coupled 
to, which is maintained through large distances (in a pure 
quasi-particle picture this distance is infinite). 

As long as both the experimental and the theoretical 
situation is not solved completely, PCS is unsuitable as 
'smoking gun'-technique to determine the symmetry of 
the HFSC order parameter. 

After finishing our work we became aware, that other 
authors have also derived the BTK-theory using an 
Hamiltonian approach to PCS. This work has been sup- 
ported by the Sonderforschungsbereich 252 Darmstadt- 
Mainz- Frankfurt, the Deutsche Forschungsgemeinschaft 
and in parts by the National Science Foundation under 
Grant No. PHY94-07194. One of us (FBA) hkes to thank 
the ITP, Santa Barbara for its hospitality. 
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